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ABSTRACT 

We present host stellar velocity dispersion measurements for a sample of 88 broad-line quasars at 0.1 < z < 1 
(46 at z > 0.6) from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. High signal- 
to-noise ratio coadded spectra (average S/N « 30 per 69 kms" 1 pixel) from SDSS-RM allowed decomposition 
of the host and quasar spectra, and measurements of the host stellar velocity dispersions and black hole (BH) 
masses using the single-epoch (SE) virial method. The large sample size and dynamic range in luminosity 
(As|oo = 10 43 2-44 7 ergs -1 ) lead to the first clear detection of a correlation between SE virial BH mass and host 
stellar velocity dispersion far beyond the local universe. However, the observed correlation is significantly 
flatter than the local relation, suggesting that there are selection biases in high-;; luminosity-threshold quasar 
samples for such studies. Our uniform sample and analysis enable an investigation of the redshift evolution 
of the M. -cr* relation relatively free of caveats by comparing different samples/analyses at disjoint redshifts. 

We do not observe evolution of the M. - er* relation in our sample up to z 1, but there is an indication that 
the relation flattens towards higher redshifts. Coupled with the increasing threshold luminosity with redshift in 
our sample, this again suggests certain selection biases are at work, and simple simulations demonstrate that 
a constant M . - cr* relation is favored to z ~ 1. Our results highlight the scientific potential of deep coadded 
spectroscopy from quasar monitoring programs, and offer a new path to probe the co-evolution of BHs and 
galaxies at earlier times. 
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The observed correlations between black hole (BH) mass 
and host galaxy properties in the local universe have moti¬ 
vated extensive theoretical work on AGN feedback and the 
co-evolutio n of supermass i ve black holes (SM BHs) and their 
hosts (e.g., iBensonl 1201 (1 iKormendv & Hoi 120131 and ref¬ 
erences therein). The M, - cr* relation between BH mass 
and bulge stellar velocity dispersion (cr*) arguably repre¬ 
sents the tightest among these relations, and has been well 
measu re d for both inactive galax ie s (e.g.. IGebhardt et al.1 
2000a! iF errarese & Merritt! 120001 iTremaine et al.1 120021 
Giilte kinet al.l 2009£_ McConnell_&Ma 20j3|) and active 


Nelson et al. 


Greene & Ho 


200 4; On ken et al. 


2006b; Dasvra etakl 


Grier et al.ll2013HWoo et al.ll2013!lBennert et al.ll2014l) in the 


galaxies (e.g . . IGebhardt et~ 2000b! FerrareseetalJ 


2004; Barth et al 


2007 


2001 


Watson et al. 


2005 


2008 


local universe. 

Given the importance of these BH scaling relations 
in our understanding of galaxy formation and evolution, 
there have been extensive observational efforts in the past 
decade both to refine the measurements of the local rela- 


Graham et al. 2001; Tremaine et al. 2002F Marconi & Hunt 

2003; Haring & Rix 2004; 

Aller & Richstone 

2007; 

Bentz et al.) 2009; 

Giiltekin et al. 

2009; Graham et al. 

2011; 

McConnell & Mai 

2013t Kormendv & Ho 2013) and refer- 


ences therein), and to probe their possibl e redshift evolution 
(e.g.. iTreu et al] 2004t iPeng et al.ll2006allbt IWoo et al.1 200 


Treu et al~T 2007: IWoo et alT~ 2008! Safvjandgr_gtak 2007[ 

Shenetal|_ 2008bLLjahnkeetaLl_ 2009^ McLeo^&Bechtold 


2009; ^ecarhetalJj20i0yMerloni_e^d. 12010! Bennertetal 

2010: ICisternas e t al.11201 1! ICanalizo et al.1 2 012! iHiner et al.. 

2012; _ Targett et al.1 l2012t Salviander & Shieldsl 12013! 

Matsuoka et al.l 12014! iBongiorno et al l 120141) . Given the 
small angular size of the sphere of influence of the BH, so 
far all measurements beyond the local universe are based 
on broad-line quasar samples, where the BH mass can be 
estimated using the properties of the broad lines, instead of 
spatially resolved gas or stellar kinematics. Most of these 
studies reported a positive evolution, advocating a scenario 
where BH growth precedes that of the bulge. However, 
statistical biases resulting from uncertainties in BH mass 
estimates and host property measurements, and from the 
selection biases in broad-line quasar samples with luminosity 
thresholds provide cautionary notes on such claims (e.g. , 
Salviander et al.l 12007! lLauer et al.l 12007! IShen et al.|[2 008t 
Shen & Kellvl 12010! iSchulze & Wisotzktll20lU l2014E ~ fShei 
201 3! ISalviander & Shieldsll201 31) . 


On the other hand, measuring the host galaxy stellar veloc¬ 
ity dispersion of quasars is notoriously difficult. It requires 
sufficient spectral quality to measure the host stellar absorp¬ 
tion features diluted by the quasar continuum. Such exercises 
based on survey-quality spectra are limited to low redshifts 
(e.g.. iGreene & Holl2006albl IShen et al.ll2008bl) . At z > 0.3, 
IWoo et al.l (120061) and IWoo et al. (120081) presented cr* mea¬ 
surements in 14 quasars at z = 0.36 and 5 quasars at z = 0.57, 
respectively, based on high S/N Keck spectroscopy. The lat¬ 
ter two studies represent the best effort to date in measuring 
cr* for higher-" quasar hosts. So far there has been no direct 
measurement of cr* in z > 0.6 quasars. 

Here we use spectroscopic data from the the Sloan Digi¬ 
tal Sky Survey Reverberation Mapping (SDSS-RM) project 
(IShen et alj|2015l) to study the relationship between BH mass 
and host stellar velocity dispersion in quasars at 0.1 < z < 1, 
a regime where BH masses can only be estimated (other 


than from BH-galaxy scaling relations^ in broad-line quasars 
with reverberation mapping (e.g.. iPetersonl 120141) or its ex¬ 
tension, the single-epoch virial BH mass estimators (e.g., 
lShenil2013l) . SDSS-RM is a major multi-object reverberation 
mapping project; it spectroscopically monitors a flux-limited 
(i < 21.7) sample of 849 quasars in a_single 7 deg 2 field with 
the SDSS-III BOSS spectrograph ( Smee et akl 120131) on the 
2.5-m wide-field SDSS telescope (Gunn et alj|2006l) . With a 
total of 32 spectroscopic epochs totaling ~ 60 hrs of exposure 
from the SDSS-RM baseline program in 2014, the coadded 
spectra provide very high S/N to allow decomposition of the 
host and quasar spectra, and measurements of BH and host 
properties. In this work we focus on cr*, and more detailed 
studies of the host stellar populations of these quasars will be 
presented elsewhere. 

The paper is organized as follows. In §|2] we describe the 
data, the host/quasar spectral decomposition, and the mea¬ 
surements of cr* and BH mass. In §0we present the results. 
We discuss the caveats of our analysis and the interpretation 
of the data in §|4]and conclude in §0 A flat ACDM cosmology 
is adopted throughout with fly = 0.3 and ho = 0.7. We use M. 
to denote the true BH mass, and to denote the virial 

BH mass estimate. By default wavelengths are restframe vac¬ 
uum wavelengths. For ease of discussion, we use “high-z” to 
refer to z > 0.3 in this paper. 

2. DATA AND SPECTRAL MEASUREMENTS 

T he spectra were taken as part of the SDSS-I II BOSS sur¬ 
vey (lEisenstein et alJl201 1! iDawson et alJl2013h from Jan to 
Jul 2014. The SDSS-RM field was observed with an upgraded 
version of the pair of multi-object fiber spectrographs used for 
SDSS I/II (ISmiee et alJ2013l) . with a fiber diameter of 2". The 
wavelength coverage of BOSS spectroscopy is 3650-10,400 
A, with a spectral resolution of R ~ 2000, and a pixel size of 
10 -4 in log 10 A (69kms _1 native SDSS pixel). The typical S/N 
per pixel in g band in a 2-hr exposure is ~ 4.5 at g pst = 21.2. 
The spectra for individual epochs were flux-calibrated using 
a custom approach based on simultaneous standard star ob¬ 
servations to an absolute accuracy of ~ 5% for point sources 
(IShen et al 120151) . 

All 32 epochs from the 2014 observations of SDSS-RM 
were coadded using an inverse-variance-weighted spline fit 
(as in the IDE-based BOSS ninclincTl idlspec2d), with a log¬ 
arithmic wavelength binning of I0 4 in log 10 A. The spectra 
are stored in vacuum wavelength, and calibrated to the he¬ 
liocentric reference with an accuracy of better than 5 kms' 1 . 
Redshifts were measured using the SDSS redshift pipeline 
dBolton et all2012i) on the coadded spectra and visually veri¬ 
fied. 

We limit the current study to the 212 z < 1.09 quasars in 
the full SDSS-RM sample for which we can decompose the 
host component from the spectrum, and measure the virial 
BH mass using the broad H j3 line, arguably the best line to 
estimate a virial BH mass (see discussion in IShenll2013L Sec 
3.1.7). All data and measurements in this paper are provided 
in a FITS catalog described in Table |T] 

2.1. Spectral decomposition 

We decompose the quasar and host-galaxy spectra us¬ 
ing the_ 2 rinci 2 alcomponentanalysis_(PCA}method adopted 
by IVanden Berk et all (120061) and IShen et aTl (l2008bl) . This 

IS Publicly available at http://www.sdss3.org/svn/repo/idlspec2d/ 











































































































































SDSS-RM: TheM.-u* Relation at High z 


3 


SDSS J141324.28+530527.0 



Figure 1 . An example of the spectral decomposition (top) and measurements of cr* (middle) and quasar properties (bottom). In the top panel, the estimated 
host fraction within the wavelength range 4160 A < A < 4210A and the redshift are given in the upper-right comer. The black line is the original spectrum, with 
spectral errors shown as the gray line. The blue and green lines are the reconstructed quasar and host galaxy spectra, respectively. The red line is the sum of the 
blue and green lines. There are residuals in the quasar emission line region due to the limited number of quasar eigenspectra used. In the middle panel, the black 
line is the decomposed galaxy spectrum (original spectrum minus reconstructed quasar spectrum) and the red line is the best-fit from vdispfit. The gray and cyan 
shaded bands are regions around quasar broad/narrow emission lines and galaxy narrow emission lines, that were masked for the fits. The best-fit cr* and the 
median S/N per pixel of the decomposed host spectrum are marked. In the bottom panel, the black line is the decomposed quasar spectrum (original spectrum 
minus reconstructed galaxy spectrum). The brown and blue lines are the model for the continuum and Fe II emission; the green line is broad H/3 and the cyan 
lines are narrow lines. The red line is the sum of all model components. The FWHM of broad H (3 is shown as the black horizontal segment within the broad line 
profile. 


method assumes that the composite (quasar+galaxy) spec¬ 
trum can be described by the combination of two indepen¬ 
dent sets of eig enspectra derived from pu re galaxy and pure 
quasar samples. IVanden Berk et all d2006l) found empirically 
that a PCA decomposition using only the first few galaxy 
and quasar eigenspectra can reasonably recover the properties 
of the host galaxy, provided there is a significant contribu¬ 
tion from the host galaxy in the composite spectrum. Using 
more eigenspectra will lead to better overall reconstruction 
of the full spectrum primarily around the broad lines. How¬ 
ever, the decomposition of the quasar and host components 
may become ambiguous because the orthogonality between 
the higher-order host and quasar eigenspectra is not guaran¬ 
teed. 

Following IVanden Berk et all (120061) and IShen et all 
(!2008bl) . we use the first five galaxy eigenspectra and the first 
ten quasar eigenspectra constructed from the SDSS (lYip et all 


12004 aiibh to decompose our coadded quasar spectra. Fig. Q] 
(top) shows an example of the spectral dec ompositi on. A de¬ 
composition is deemed to be successful bv IVanden Berk et alJ 
(120061) if the host fraction fa within the wavelength range 
4160A < A < 4210 A is greater than 0.1. This threshold 
of course depends on the quality of the spectrum. For our 
high S/N spectra, we set a threshold of fa = 0.05 for a 
successful decomposition, with additional criteria to test 
the reliability of the successive spectral measurements on 
the decomposed components. As in IShen et al . (120081 4). we 
obtain the pure-quasar (pure-host) spectrum by subtracting 
the reconstructed galaxy (quasar) spectrum from the original 
spectra, and use the original spectral errors for the decom¬ 
posed spectra. We refer to these spectra as the (decomposed) 
quasar and host spectra. These spectra are used in our fitting 
procedures instead of the directly reconstructed spectra, 
because the former allows use of the associated spectral 
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Figure 2. Top: Distribution of quasars in the luminosity-redshift plane. The 
continuum luminosity is estimated at restframe 5100 A from the total spec¬ 
trum. The full sample from SDSS-RM is shown in black open circles, and 
quasars with good cr* measurements are shown in r ed circ les. The gray 
dots are the low-z SDSS quasar sample inIShen et al. H2008bl) , and the cyan 
open triangles are the samples at z = 0.36 4Woo et all12006 ) and at z = 0.57 
(Woo et al. 2008), as we will discuss in a Bottom: Redshift histograms for 
the parent sample (black) and the sample with good cr* measurements (red). 

errors, and the regions with large residuals due to imperfect 
decomposition can be easily masked during the fitting. In 
addition, the reconstructed spectrum is a “model” rather than 
real data, and there could be systematic differences between 
the reconstructed spectrum and the real spectrum due to the 
limited number of eigenspectra used. 

The spectral decomposition removes broad-band quasar 
features such as the power-law continuum and Fe II emission, 
as well as narrow-band features including quasar broad and 
narrow emission lines from the galaxy spectra. While the 
removal of the broad-band quasar feature is reasonably ac¬ 
curate, the removal of the quasar emission lines often shows 
significant residual, due to the small number of quasar eigen¬ 
spectra used in the decomposition (which are insufficient to 
account for the full variance seen in quasar emission lines). 
We could add higher-order quasar eigenspectra to better re¬ 
construct the quasar broad and narrow emission lines. But 
doing so would potentially over-fit the quasar component in 
the stellar absorption line regions and confuse the decompo¬ 
sition, which will impact the cr* measurement. Therefore we 
opt to use a small number of quasar eigenspectra to remove 
most of the broad-band quasar features, and then mask out 
these emission-line regions when measuring cr* with the de¬ 
composed host spectrum. 

2.2. Stellar velocity dispersion measurements 

The decomposed galaxy spectrum can be used to measure 
properties of the stellar population in the quasar host. In this 
work we focus on the stellar velocity dispersion, as the mea¬ 
surement of this quantity is relatively less affected by the im¬ 
perfections in the spectral decomposition compared to other 
properties of the stellar population. We defer a more detailed, 
full analysis of the host stellar populations to future work. 


We use the vdispfit routine, publicly available in the idl- 
spec2d product, to measure cr* from the decomposed galaxy 
spectra. This routine is used to measure stellar velocity dis¬ 
persions of SDSS galaxies in the standard SDSS spectro¬ 
scopic products, vdispfit is a template fitting algorithm that 
finds the best-fit cr* by minimizing the \ 2 in pixel space. 
The templates consist of the first four PCA eigenspectra con¬ 
structed from the echelle stellar spectra in the Elodie database 
(iMoultaka et al.ll2004l) and cover restframe wavelengths from 
4125 A to 6796 A. 

We restrict the fit to the wavelength range of 4125-5350 A, 
which includes copious stellar absorption features such as 
the G band (4304 A), the Mg lb AA 5167,5173,5184 triplet, 
and Fe (5270 A), but excludes the Ca H+K AA 3969,3934 re¬ 
gion that is inherently unreliable for stellar velocity disper¬ 
sion measurements (e.g.. iGreene & Holl2006a f. In addition, 
we mask narrow emission lines due to both galaxy and quasar 
excitation with a velocity window of ±400kms' 1 . The spec¬ 
tral region around the quasar broad H/3 line (4760-5020 A), 
where the PCA decomposition often shows significant resid¬ 
uals, is also masked for the fits. A small wavelength shift 
(within ±500kms -1 ) is allowed in the fit to account for the 
uncertainty in the systemic redshift estimate. A fifth-order 
polynomial is added to account for the broad-band continuum 
shape. Fig.|T|(middle) shows an example of the vdispfit fits to 
measure cr*. The reported cr* by vdispfit is already corrected 
for the native SDSS instrument broadening (~ 65kms -1 ). 

As a sanity check, we also use the direct-pixel-fitting code 
Penalized Pixel-Fitting (pPXF) of ICappellari & Emselleml 
(2004) to measure the stellar velocity dispersion. We use the 
stellar library from IValdes etal . (120041) with a spectral res¬ 
olution of 25 kms -1 and correct the templates to the SDSS 
spectral resolution before fitting. Only the first two moments 
(radial velocity and velocity dispersion) are fitted. In addition 
to the radial velocity and velocity dispersion, the code solves 
for an eighth order multiplicative polynomial that accounts 
for flux calibration errors. Errors in the dispersions are de¬ 
rived by generating 100 mock spectra as a sum of the original 
spectrum and a Gaussian random realization of the spectral 
noise. Only the spectral range of 4000-5350A is used in the 
fit. We mask out regions around quasar narrow but not broad 
emission lines. The differences in the code, spectral fitting 
range, and masks between the pPXF and vdispfit approaches 
are designed to assess the systematic uncertainty in cr* mea¬ 
surements. 

The success of the cr* measurements depends on a number 
of factors such as the subtraction of the quasar spectrum, the 
S/N of the resulting host galaxy spectrum, and the strength of 
the stellar absorption features. We define the following set of 
criteria for a successful cr* measurement: (1) the host fraction 
f H is greater than 0.05; (2) the median S/N per pixel of the de¬ 
composed host galaxy spectrum over restframe 4125-5350 A 
is greater than 4; (3) cr* is measured at > 3cr statistical confi¬ 
dence level; and (4) visual inspection shows that the fit is not 
severely affected by poor PCA decomposition (i.e., significant 
residuals exist in the fitting region). The first three criteria re¬ 
move 104 objects and the last criterion removes 20 additional 
objects, all with no prior knowledge of BH mass. 

When compare the resulting 88 cr* values measured with 
both vdispfit and pPXF, there were 18 objects for which the 
difference is larger than 2 cr. However, 16 of these 18 ob¬ 
jects are fits with small statistical errors in high S/N spec- 
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tra, and the fractional difference between the two methods 
is < 30%. Hence the difference between the two methods 
can be ascribed to the underestimation of the systematic un¬ 
certainties in the cr* measurements (i.e., due to differences in 
stellar templates, fitting ranges and emission line masks). We 
therefore retain the cr* values from the fiducial vdispfit method 
and indicate by the SIGMA_WARNING flag in Table Q] that 
their measurement errors may underestimate the true uncer¬ 
tainty. For the remaining two objects, manual inspection sug¬ 
gests that the pPXF result is more robust (based on a smaller 
reduced x 2 ) and we decided to replace their vdispfit cr* mea¬ 
surements with those from the pPXF method. Our final sam¬ 
ple of quasars with good cr* measurements consists of the 88 
objects (which we call the “good sample”). Fig.[2]shows their 
distribution in luminosity-redshift space as compared to the 
parent sample. Although objects at high z are preferentially 
lost due to our fn and S/N cuts, there is otherwise no appar¬ 
ent bias in the good sample compared to the parent sample in 
terms of object luminosity. In particular, the success rate of 
cr* measurement is quite high at z < 0.75 (~ 70%), ensuring 
that our good sample is representative of the general quasar 
population at these redshifts. 

Although we have intentionally avoided the Ca H+K region 
in our cr* fits, these features are the strongest absorption fea¬ 
tures in most of our objects, and therefore easier to measure 
than those at longer wavelengths. For the sake of complete¬ 
ness, we also measured cr* in the restframe 3780-4000 A re¬ 
gion using pPXF (the vdispfit templates do not extend to this 
regime), and compile the measurements in Table Q] There 
is broad agreement between these Ca H+K dispersions and 
our fiducial dispersions, but the scatter is large. We do not 
use these Ca H+K dispersions in our following analysis, and 
merely provide them in Table Q] for future reference. 


2.3. Single-epoch virial mass estimates 

To apply single-epoch virial BH mass estimators for our 
quasars, we follow earlier work (e.g., iShen et alJ200M[2()nh 
to fit the quasar spectrum with continuum emission, Fe II 
emission, and broad and narrow emission lines around the 
H j3 region. Fig. Q] (bottom) shows an example of the fits to 
the quasar spectrum. The continuum luminosity (i.e., L 5100 = 
AL A 5100 j j 1 ) and broad line width are measured from the best- 
fit model. We use the Monte Carlo approach described in 
IShen et aT] (1201 ll) to estimate measurement errors for these 
quantities. We choose the FWHM to characterize the width of 
the broad H/i line instead of using the line dispersion (j\ mc for 
the reasons discussed in lShenl ( 2013 1 and references therein). 
However, we will experiment with different virial BH mass 
recipes that account for the difference between using FWHM 
and cr \[ ne . 

The standard SE virial BH mass estimator based on H/3 is: 


log 


( Tf B H,vii- 

V 


= a +/? log 


-A 100 


10 44 erg: 


,-1 


+ clog 


/ FWHM \ 
\ kms ' 1 ) 


(1) 

where coefficients a, b and c are determined empirically us- 
ing the local RM AGN sample in different recipes (e.g., 
Vestergaard & Peterson! [2006 ; IWang et~aki [20091: IFeng et afl 
20141) . The luminosity range in the local RM AGN sample 
encompasses most of our quasars. 

Our fiducial mass recipe is that of IVestergaard & Petersonl 
(120061 eqn. 5; hereafter VP06) with a = 0.91, b = 0.5 and c = 2. 
In VP06, the continuum luminosity L 5100 was not corrected 



o* [km s' 1 ] 

Figure 3. The relation between cr* and SE virial BH mass for the sample 
with good velocity dispersion measurements. The black points show all 88 
objects, and the red points show the 46 objects at z > 0.6. Error bars are mea¬ 
surement errors only, and the systematic uncertainty in the virial BH masses 
is ~ 0.5 dex < fShcnll2()T3l >. Spearman's test results are shown: a significant 
correlation is detected both for the full sample and for objects aL z. > 0.6. The 
black solid line is the median relation derived using the Bayesian linear re¬ 
gression method developed by IKellvl (2007 ). and the gray band indicates the 
2a confidence range of the correlation. The vertical scatter around the best-fit 
relation, after correcting only for measurement errors, is 0.4 dex. The dashed 
line is the relation for local inactive galaxies from IKormendv & Hd (2013 . 
eqn. 7). The slope of the MgH.fir - relation in our quasars is significantly 
shallower than that of the local relation, as discussed in the text. 

for host starlight contamination, hence we use the total 5100 A 
luminosity in this red pc FI We use the VP06 virial BH masses 
in our following analysis. 

We also consider an alternative recipe in IFeng et al.l d2014l 
eqn. 14; hereafter F14) with a = 3.602, b = 0.504 and c = 1.2. 
The F14 recipe requires the host-corrected quasar Tsioo.qso, 
which is measured from the decomposed quasar spectrum. 
This recipe relaxes the assumption that BH mass is propor¬ 
tional to the square of FWHM, and uses the best-fit slope on 
FWHM from the calibration against^the local sample of AGN 
with RM masses (also see, e.g JWang et al.l2009l) . Since these 
RM masses were computed using the line dispersion crime in 
the rms spectrum rather than FWHM, this recipe effectively 
accou nts for the nonlinear relation betwee n FWHM and crii ne 
(e.g JCollin et alJl2006l : IShen & Holl2014l) . We use this alter¬ 
native set of BH mass estimates to demonstrate the ambigui¬ 
ties in measuring BH masses for high-; quasars with the SE 
methods, as discussed in ^4] 

The measurement errors in logM B H,vir are propagated from 
the measurement errors in FWHM and logZ^oo, but the dom¬ 
inant uncertainty in logMBH,vir is the systematic uncertainty 
of single- epoch virial masses, which can be as large as ~ 0.5 
dex (e.g.. lShenll2013l) . We successfully derived Mbh.vL esti¬ 
mates for all quasars in our parent sample except for a few 
cases where restframe 5100 A is not covered by the spectrum. 

15 We also tried using the quasar-only 5100 A luminosity in the VP06 
recipe and found similar conclusions. The systematic uncertainties in the 
virial BH mass estimates (including the zero-point uncertainty) are much 
larger than this detail (e.g.. lShenl20l3i ). 
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Table 1 

Catalog format 


Column 

Format 

Description 

RMID 

LONG 

Index in the SDSS-RM master catalog [0-999] 

SDSS_NAME 

STRING 

SDSS designation based on RA/DEC (hhmmss.sszbddmmss.s) 

RA 

DOUBLE 

J2000 Right ascension [degree] 

DEC 

DOUBLE 

J2000 Declination [degree] 

Z 

FLOAT 

Redshift 

WAVE 

DOUBLE[4649] 

Observed vacuum wavelength [A] 

FLUX 

DOUBLE[4649] 

Total flux density [10 -17 ergs _1 cm -2_1 ] 

ERR 

DOUBLE[4649] 

Flux density error [10 _l7 erg s -1 cnT 2-1 ] 

FLUX_GALt 

DOUBLE[4649] 

Decomposed host galaxy spectrum [10 _17 ergs _1 cm _2 A ] 

FLUX_QSO t 

DOUBLE[4649] 

Decomposed quasar spectrum [10 _17 ergs _1 cm _2 A ] 

MEDSN_TOT 

FLOAT 

Median S/N per pixel for the full spectrum 

MEDSN_GAL 

FLOAT 

Median S/N per pixel for decomposed galaxy spectrum over restframe 4125-5350 A 

F_H 

FLOAT 

Host fraction within restframe 4160^-210 A; 0 if decomposition failed 

F_H_5100 

FLOAT 

Host fraction at restframe 5100 A 

SIGMA 

DOUBLE 

Fiducial stellar velocity dispersion [km/s]; 0 if not available 

SIGMA_ERR 

DOUBLE 

Measurement error in SIGMA; —1 if not available 

SIGMA_ERR_WARNING 

LONG 

1 if SIGMA_ERR may underestimate the systematic uncertainty 

SIGMA_HK 

DOUBLE 

Stellar velocity dispersion based on Ca H+K [km/s]; 0 if not available 

SIGMA_HK_ERR 

DOUBLE 

Measurement error in SIGMA_HK; —1 if not available 

LOGL5100_TOT 

DOUBLE 

Total continuum luminosity at restframe 5100 A [erg/s]; 0 if not covered 

LOGL5100_TOT_ERR 

DOUBLE 

Measurement error in LOGL5100_TOT 

LOGL5100_QSO 

DOUBLE 

Quasar continuum luminosity at restframe 5100 A [erg/s]; 0 if not covered 

LOGL5100_QSO_ERR 

DOUBLE 

Measurement error in LOGL5100_QSO 

FWHM_HB 

DOUBLE 

FWHM of broad H/3 [km/s] 

FWHM_HB_ERR 

DOUBLE 

Measurement error in FWHM_HB 

LOGMBH_VP06 

DOUBLE 

SE virial mass estimate based on the VP06 H f3 recipe [Mq] 

LOGMBH VP06 ERR 

DOUBLE 

Measurement error in LOGMBHJVP06 


NOTE. — Format of the fits catalog containing the data and measurements of the 212 quasars used in this paper, ordered by RMID. The 88 
quasars with good cr» measurements are identified with SIGMA > 0. The SIGMA_ERR_WARNING flag is set if the difference between 
the fiducial cr* measurement and that from pPXF is larger than 2a. The errors reported in SIGMA_ERR may be an underestimation in these 
cases. The <r* measurements based on Ca H+K features (SIGMA_HK and SIGMA_HK_ERR) are not suitable for studying the M. -cr* 
relation, and are only included here for future reference. The full catalog is available in the electronic version of the paper. 

^The decomposed galaxy (quasar) spectra are obtained by subtracting the reconstructed quasar (galaxy) model from the full spectrum, and 
onl y cov er the rest-frame 3450-7498 A range, i.e., the overlapping wavelength range in the quasar and galaxy PCA eigenspectra from 
I Yip et alj 120043) and lYmetaLl d2004bl) . 


None of these sources were included in the good sample. 

The VP06 SE mass recipe uses a constant average virial 
coefficient that normalizes the virial BH masses of nearby re¬ 
verberation mapping AGN to follow the local M. - cr* rela¬ 
tion in inactive galaxies (e.g. JOnken et al.ll2004l) . An implicit 
assumption in our investigations is that this average virial co¬ 
efficient does not evolve with redshift, which is a reasonable 
assumption since this is largely a geometrical factor to ac¬ 
count for the structure of the quasar broad-line region (e.g., 
IShenll20l3h . 


3. RESULTS 

We show our cr* and log/V/ bh.vh measurements in Fig. [3 
The full sample is shown in black circles with measurement 
errors only. The Spearman test finds a strong probability of a 
correlation, as shown by the values presented in Fig. [3] For 
the first time, we were able to detect directly a correlation 
between cr* and virial BH mass estimates at (z) = 0.6, while 
earlier studies did not have the statistics and dynamic range to 
clearly detect a correlation at z > 0.3 (e.v.- lWoo et all120061. 
2008). A moderately strong correlation is still present if we 
restrict the analysis to the subset of 46 z > 0.6 objects (red 
circles). The correlation for the high-z subset is similar to the 
relation for the full sample, suggesting negligible evolution, 
as further demonstrated below. 


Table 2 

Linear Regression Results 


Sample 

(z> 

a 

p 

Intrinsic scatter 

full 

0.60 

q q'7'7+0.067 

O-3//-0.066 

1 535+ 0304 
-0.303 

0 406+ 0035 

u.+uo_q 03i 

z > 0.6 

0.76 

o onc+0.088 

o. 0 88 

1 OR 1+0-511 
i.Uol_o 492 

0 440 +0 059 
U.HH-U_oo48 

0.12 <z< 0.42 

0.26 

8-3247!!;! ™ 

1 6QS+ 0 - 614 
Loyj - 0 . 604 

0 4 io +0 082 

U.H-IU_0 063 

0.42 < z < 0.60 

0.53 

O 0 79+0.137 
°‘ J ' -0.138 

1 SQ7+0-677 

1 -0.660 

0 414+ 0 082 
-0.062 

0.61 < z < 0.76 

0.70 

0 OOO+0.233 
O.3O5_ 0 235 

| iqq+1.231 

LLyy -l .216 

0 536 +0 109 
u,JJU -0.0 83 

0.76 <z< 1.00 

0.84 

8.36478:811 

0 960 +0 - 598 

u.you_Q 588 

H 981+0.091 

U.95 l_o 068 


NOTE. — Linear regression of the relation log(MsH,vir/M 0 ) = 
a + ^ilog(cr,/200kms _I ) on the data, using the Bayesian ap¬ 
proach described in iKellvl 12007 ). The intrinsic scatter refers 
to the scatter in logA/stt around the best-fit relation, after only 
accounting for measurement errors. Systematic uncertainties in 
logMBH.vir are included in this intrinsic scatter. The lcr errors 
are estimated from the 16% and 84% percentiles of the posterior 
distribution of the parameter from the Bayesian fits. 

We use the Bayesian linear regression method in Kellvl 
(120071) to fit for logMeH. v i r as a function of log cr*: 


log 


Mbh,\ 

Mo 


= a + /31og( 


200 km s 


( 2 ) 
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Figure 4. Redshift evolution of the Mbh.vu- _ <?* relation. The full sample (88 objects) is divided into four equal-number subsamples by redshift. For each 
subsample, we perform the same Bayesian linear regression on the data. In each redshift bin, the red line is the median relation from the Bayesian fit; the cyan 
and light blue shaded regions are the lcr and 2 a confidence ranges of the fit. For reference, the median relation for the full sample is shown as the black solid 
line, and the local relation in IKormendv & Hoi 42013D is shown as the black dashed line. We do not detect significant differences in the MBH,vir — 0* relation from 
the lowest redshift bin to the highest redshift bin in o ur sample. Howev er, there is an indication that the relation becomes flatter towards higher redshift, which is 
further discussed in 0 We also show the samples in I Woo et all d20061 z = 0.36) and IWoo et all 42008L z = 0.57) in red open triangles. The objects in Woo et al. 
are more luminous than our quasars in the relevant redshift bins (see Fig. [2)- Since Woo et al. did not report FWHM measurements, we took BH mass estimates 
directly from their papers. 


log(M BH vil /M©)=a+piog(G*/200 km s' 1 ) 



0.2 0.4 0.6 0.8 

Redshift 


Figure 5. Redshift evolution of the best-fit parameters in the MsH.vir — cr* 
relation defined by Eqn. m and tabulated in Table [2] The black filled circles 
are for the full sample, the red open circles are for the z > 0.6 subsample, and 
the cyan open squares are for the four redshift bins. 

where a and /3 are the intercept (at cr* = 200 kms -1 ) and slope 
of the regression, and we use e to denote the intrinsic scat¬ 
ter in logM BH around the best fit. This Bayesian regression 
method properly accounts for the intrinsic scatter, thus the 
best-fit slope is not biased due to a few points with small mea¬ 
surement errors (e.g., see discussion in lTremaine et al.lf2002l) . 
The results are summarized in Table [2] where the significantly 
non-zero (at 5cr) slope 3 for the full sample again suggests 
a strong correlation. We only consider the measurement er¬ 


rors of logMBH,vir in the fit, and thus the intrinsic scatter e in¬ 
cludes the contribution from systematic errors in logM BH:V ir- 
We also tried adding 0.5 dex systematic errors on logM B H.vii , 
or using alternative regression methods (e.g.. lTremaine et all 
120021) . but did not find significantly different results for the 
intercept and slope. 

In Fig. [3 the solid line represents the relation with the me¬ 
dian intercept a = 8.38 and slope [3= 1.54 in the posterior dis¬ 
tribution from the Bayesian fit, and the gray band indicates the 
2a confidence range of the correlation. The intrinsic scatter 
is estimated to be 0.4 dex. Although there is a strong correla¬ 
tion between M B H, v ir and cr* in our sample, the best-fit slope 
is significantly shallow er than that of the M, - c r* relation for 
local inactive galaxies (Korm endv & Hoi 120131 eqn. 7, with 
j3 = 4.38), which is shown as the dashed line. We discuss this 
discrepancy in detail in ®] On the other hand, the average BH 
mass at cr* = 200 kms -1 is ~ 2.5 x 10 8 Mq, consistent with 
the local relation. 

The large dynamic range in luminosity and redshift of our 
sample, as well as the homogeneous analysis, make it possible 
to investigate potential redshift evolution of the M KHvu - a * 
relation within our sample. We divide the full 88 objects into 
four equal-number subsamples by redshift, and show the re¬ 
sults in Fig. [4] We perform linear regression with the Bayesian 
method to the data in each subsample, and show the median 
relation (red lines) and the lcr (2cr) confidence range in cyan 
(light blue) shaded bands in Fig. [4] A correlation is inferred 
(albeit at low statistical significance) even in the highest red¬ 
shift bin with (z) = 0.84. The regression results are summa¬ 
rized in Table [2] and plotted in Fig. [5] 

The median relation for the full sample is shown in black 
lines, and is consistent within lcr with the relations at differ¬ 
ent redshifts. Therefore we do not detect any redshift evolu¬ 
tion in the M B H.vir ~ cr* relation. In particular, the mean BH 
mass at cr* = 200 kms -1 is remarkably stable, and js consis- 
tent w ith the value from the local relation in IKormendv & Hoi 
(120131) . There is, however, an indication that the relation flat¬ 
tens towards higher redshift (see Tabled, as further discussed 
in®] 

4. DISCUSSION 

The results presented above should be interpreted with 
great caution. Below we discuss a list of potential caveats 
that went into the above analyses and interpretation, and our 
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l°g L 5ioo,tot 43.0 43.5 44.0 44.5 



Figure 6. Comparison between our results and earlier results. All data points 
are color-coded by the total continuum lumin osity L smo- Measurement errors 
have been suppressed for clarity. The IShen et ahl (200BH) sample includes 
quasars at lower redshift, and their luminosities are substantially lower than 
our quasars, as indicated in Fig. [2] The lShen et all 12008N) sample shows sys¬ 
tematic offset towards lower BH masses, which is an indication that our high- 
Z sample is mis sing low-mass BH s at fixed i t* due to the luminosity threshold. 
The samples in Woo et al. (2006) and Woo et al. 12008) are shown in open tri¬ 
angles. They overlap with our quasars with comparable luminosities, but are 
limited to a narrow dynamic range in cr* and BH mass. Since Woo et al. did 
not report FWHM measurements, we took BH mass estimates directly from 
their papers. The black line and gray shaded band are the Bayesian linear 
regression results as in Fig. [3 and the dashed line denotes the local relation 
from IKormendv & Hoi 12013k 

resolutions. 

Spectral decomposition 

Measuring the host stellar velocity dispersion in quasars 
is always challenging because the stellar absorption features 
are diluted by the significant contribution from the quasar. 
The PCA decomposition technique attempts to subtract the 
quasar continuum and broad line emission (including the 
optical Fe II complex that extends over a broad spectral 
range). However, it is difficult to ensure a clean subtraction 
of the quasar contribution, and this may affect the dispersion 
measurements. To mitigate the residual impact of quasar 
emission, we masked prominent broad and narrow emission 
line regions when fitting the decomposed galaxy spectrum for 
er*; while broad-band quasar residuals are accounted for with 
low-order polynomials in the er* fit. We believe in general 
this approach should work well. 

er* measurements 

The main hurdle in our er* measurement is the S/N of the 
stellar absorption features in the decomposed galaxy spec¬ 
trum. We have set a threshold of median S /N > 4 per pixel 
to ensure a reasonable er* measurement (the median S/N per 
pixel in the galaxy spectrum is ~ 15 for our good sample). 


We have also tested the sensitivity of the er* measurement to 
changes in the stellar templates, fitting code, spectral fitting 
range and masks, and found our fiducial er* measurements are 
robust against these details. The ultimate limitation on er* 
measurements in quasars likely lies in the systematics from 
spectral decomposition in extremely quasar-dominated spec¬ 
tra, spectral calibration and template matching - our measure¬ 
ments are far from reaching that limit. For example, it is still 
difficult to measure cr* in the most luminous PG quasars at 
low z (which are highly biased relative to the general quasar 
population in terms of luminosity and radio-loudness) with 
a huge quasar-host contrast - S/N is not the limiting factor 
thereQ 

The completeness in er* measurements for quasar samples 
poses another concern. If a significant fraction of objects with 
high quasar-to-host contrast are preferentially lost from the 
sample due to the failure of er* measurements, there is a po¬ 
tential selection bias in the resulting M KH ylr -a t relation, such 
that the average cr* at fixed BH mass may be biased high. 
Fortunately for our sample, the completeness in cr* measure¬ 
ments is high (« 70%) at z < 0.75, meaning that the potential 
bias due to sample incompleteness is not a severe issue for 
the bulk of our sample. However, it may become important 
for the highest redshift bin shown in Fig.|4] where the fraction 
of quasars with measured er* is low mainly due to small host 
fractions. There is no easy way to quantify this effect - deeper 
spectroscopy is required to increase the chance of measuring 
er* in these objects. 

The next question is what er* we are actually measuring. 
SDSS-BOSS spectra collect light enclosed in a 2" diameter 
circle on the sky, which translates to a physical radius of 
3.3— 7.8kpc at 0.2 < z < 0.9. Therefore our er* is measured 
from the integrated light from the bulge and a potential disk 
component at high z. We are ignorant about the morpholog¬ 
ical types of our quasars and their redshift dependence, so 
an object-by-object correction to the bulge er* is difficult. 
We have tried the empirical correction for aperture effects 
in IShen et ahl (l2008bl) . and did not find significant changes 
in our results. While we caution that some of the scatter in 
our results may arise from the failure to resolve the host light 
into disk and bulge components, the estimated bias in er* due 
to the inclusion of a rotating disk is < 16% based on local 
AGN samples (e.g. jBennert et akll2014l) . Thus this effect is 
negligible given the measurement uncertainties and intrinsic 
scatter in Fig. [3] 

BH mass estimates 

A larger concern is in the BH mass estimates based on the 
SE methods. As discussed at length in IShenl (1201 3l) . there 
are currently large systematic uncertainties in the SE virial 
masses for high-z quasars at the ~ 0.5 dex level (both in scat¬ 
ter and in zeropoint), and systematic differences between dif¬ 
ferent recipes. 

To investigate the difference due to alternative SE mass 
recipes, we have repeated the analysis based on the F14 mass 
recipe for H/3. We still find a significant correlation between 
A/BH,vir and cr* in different redshift bins, but the slope of the 
correlation is even shallower than that based on the VP06 

20 A promising approach for measuring cr* in these most luminous quasars 
is to use near-infrared stellar absoiption features assisted with adaptive optics 
to minimize the impact of the quas ar continuum (e.g., Dasyra et all 120071 : 
Watson et al. 2008; Grier et al. 2013). But even in this case, the measurement 
of cr* can be extremely difficult in the most luminous quasars. 
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masses. This exercise highlights the ambiguity in estimat¬ 
ing BH masses using the SE method, and the difficulty in 
comparing different studies using different BH mass recipes. 
Nevertheless, the sm all samples at z = 0.36 and z = 0.57 in 
IWoo et all ( 2006 ) and lWoo et al. ( 2008 ). when matched in lu¬ 
minosity and redshift to our quasars, fall consistently within 
our data points in Fig. [4] 

Since all SE mass recipes explicitly depend on luminosity, 
statistical biases will arise from samples with a threshold 
luminosity (as we further discuss below). 

Statistical biases 


There are several statistical biases due to quasar-based 
sample selection and BH mass uncertainties that have 
been previously discu s sed (e.g.. lAdelberger & Steidell 12005 : 
Salviander et alJ 120071 lLauer et ah T2007 : Shen et afT l2008a: 


Shen & Kellvl 120 lof ISchulze & Wisotzki 1201 ll) . Below we 

discuss these effects, and provide qualitative assessments of 
their implications on the interpretation of Fig. [3] and Fig. [4] 
lLauer et al.l d2007l) po i nted out that (also see 
lAdelberger & Steidell 120051: ISalviander et al.l 120071) . be¬ 
cause there is intrinsic scatter in the M m - cr* relation, samples 
selected based on a threshold in quasar luminosity will prefer¬ 
entially select overmassive BHs with respect to galaxy stellar 
velocity dispersion, assuming that more lumin ous quasars 
represent more massive BHs on average. IShen & Kellvl 
(120101) further sugges ted that (based on the earlier suggestion 
in Shen et all l2008ah . because SE virial mass estimates 
explicitly depend on luminosity and there is non-negligible 
random scatter between instantaneous luminosity and line 
width (i.e., the origin of the uncertainties in SE virial masses), 
quasar samples with a threshold luminosity are, on average, 
biased high in their SE virial mass estimatesFI These two sta¬ 
tistical biases will work independently to produce a BH mass 
offset relative to cr* in high-z, hence high-luminosity, quasar 
samples; this BH mass offset is more prominent at lower cr* 
(hence flattening the slope in M.-cr*), where under-massive 
BHs (relative to cr*) are more easily lost due to the sample 
luminosity threshold. In addition, both effects become worse 
at higher threshold luminosities, which qualitatively explains 
the flattening in the slope of the log A7 BH v i r - cr* relation with 
redshift seen in Fig. [4] A corollary of this result is that if 
we probe sufficiently faint quasar luminosities at high z, we 
should recover the population of undermassive BHs at fixed 


O *. 

To strengthen these arguments, we perform a test using 
the low-r quasar sample in iShen et aTI (2008b|). Under the 
premise that there is no redshift evolution in the intrinsic 
M. — cr* relation, we plot the sample in IShen et aD (l2008bl 
hereafter the Shen08 sample), and compare with our data in 
Fig. [6] The Shen08 sample is a low-z SDSS quasar sam¬ 
ple, analyzed using the same PCA method to decompose the 
spectrum, and measure host stellar velocity dispersion, quasar 
continuum luminosity and broad H/3 FWHM. To be consis¬ 
tent with our analysis, we use the total continuum luminos¬ 
ity Lgioo from their compiled measurements, and scale their 
FWHM values upwards by a constant 0.05 dex to match our 
own measurements for the objects in common with the cat¬ 
alog of IShen et al.l (1201 ll) . We then use the VP06 recipe to 


21 A detailed derivation of this virial BH mass bias is presented in Sec 3.3 
of ISheni 120 131) . which is referred to as the luminosity-dependent bias in SE 
virial BH masses. 



too 100 100 

a, [km s ] 


Figure 7. Illustration of the statistical biases discussed in 0 We simulate 
a mock distribution of true B H masses using the observed local stellar ve¬ 
locity distribution function (Bernardi et al. 2010), convolved with the local 
Mm — cr* relation ( Kormendy & Ho 2013]) with a scatter of 0.44 dex in Mm 
(e.g. JGultekin et aD120091 We then assign quasar luminosities assuming a 
lognormal Eddington ratio (A) distribution with a mean of (log A) = — 1 and a 
scatter of 0.3 dex. We th en generate F WHMs and virial BH masses following 
the approach detailed in lShenl 420131) . The distribution of FWHM is lognor¬ 
mal with a scatter of 0.18 dex, and a mean depending on the true BH mass 
and the mean luminosity at fixed true mass. By design, the simulated virial 
BH masses have an uncertainty of 0.39 dex around their true masses, roughly 
consistent with current estimates (e.g.. lShenl2013D . We then impose different 
threshold luminosities to select a random set of 500 quasars from the mock 
sample, shown in open circles. The top panels show the distribution in terms 
of true BH mass, and the bottom panels show the corresponding luminosity- 
threshold mock samples in terms of virial BH ma sses. The sol id lines repre¬ 
sent the best-fit linear regression results using the lKellvl 420071) method. Th e 
black dashed lines denote the local relation from I Kormendy & Hd 420131) . 
As the luminosity threshold increases, we preferentially select over-massive 
BHs at the low-cr* end, which causes a flattening in the regression slope. The 
top panels demonstrate the lLauer et aD 420071) bias, wh ile the bottom pan- 
els demonstra t e the combined effect of the[Laueret aD 120071) bias and the 
IShen & Kellvl 1201(1) bias. 

compute virial BH masses for their sample. The Shen08 sam¬ 
ple is shown in Fig. [6] with small circles. This low-z sample 
has significantly lower luminosities than our quasars (see Fig. 
0, and as expected, they are systematically offset from our 
quasars at fixed cr*. In fact, a systematic trend of BH mass 
offset as luminosity increases is present in the Shen08 sample 
as well. These results are consistent with the above argument 
that high-z luminosity-threshold samples suffer from statisti¬ 
cal biases in their average BH mass relative to the expectation 
from the mean M m — cr* relation. 

To further illustrate the effects of these two statistical bi¬ 
ases, we perform simple simulations of luminosity-threshold 
quasar samples. We start by simulating a distribution 
of quasars with known velocity dispersions and true BH 
masses using the local stellar velocity dispersion function 
dBernardi et al. 2010) convolved with the local M. — cr* rela¬ 
tion dKormendv & Holl2013l) . We assume an intrinsic scatter 
of 0.44 dex in M,_ at fixed cr* for all morphological types of 
galaxies (e.g. JGiiltekin et al,l2009l) . With the true BH masses, 
we generate quasar bolometric luminosity L ho | using a log¬ 
normal distribution of Eddington ratio, A = Lboi/TEdd, where 
Z^Edd = 1 -26 x 10 38 (M./Mq) erg s -1 is the Eddington luminos¬ 
ity. We choose a mean Eddi ngton ratio (log A) = -1 and a 
scatte r ctl = 0.3 dex (e.g., IShen et all 12008a: S hen & Kellvl 
120121). W e assum e Lboi = IOL 5100 to get £ 5100 - We then fol- 
low lShenl (1201 3l) to generate FWHM and virial BH masses 
using the VP06 recipe. The FWHMs are generated assum¬ 
ing a lognormal distribution, with the mean depending on 
true BH mass and mean luminosity at fixed true mass, and 
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a scatter of ofwhm = 0.18 dex. This approach ensures that 
the virial BH masses are on average unbiased (when no lumi¬ 
nosity constraints are imposed), and that the FWHM distribu¬ 
tion of quasars is consistent with observation (e.g. JShen et all 
I2008al) . We assume no correlation between the luminosity 
dispersion <7l and FWHM dispersion ctfwhm at fixed true 
BH mass. Therefore the uncertainty of the virial masses is 
er v ir = \J (0.5ctl) 2 + (2crFWHM)i = 0.39 dex, roughly consistent 
with expectation (e.g., IShenll2013l) . In addition, this approach 
yields a maxim um luminosity-dependent bias in the virial 
BH m asses (e.g. JShen et al.ll2008atlShen & Kellvf2010HShenl 
120131) . Again, we refer the reader to iShenl (120 1 31 Sec 3.3) for 
a more elaborate derivation of this virial BH mass bias. 

The results are shown in Fig. |7] where we select a random 
set of 500 quasars from the simulated sample by imposing 
different luminosity thresholds. We have shuffled the cr* val¬ 
ues by 0.1 dex to mimic measurement errors. The top panels 
show the simulated distributions in terms of true BH masses, 
which demonstrate the effect of the lEauer et al.l (120071) bias. 
The bottom panels show the simulated distributions in terms 
of virial BH mas ses, which demonstrate the combined effect 
of the lEauer et ail ( 2007 ) bias and the virial BH mass bias. 
In both cases, the luminosity-threshold samples preferentially 
select overly-massive BHs at the low-cr* end, and the best- 
fit regression slope (the solid lines) is significantly shallower 
than the intrinsic slope (the dashed lines). This simple simula¬ 
tion qualitatively explains the observed trends in our data, and 
p rovides strong evidence that statistical biases are at work. 

iSchulze & WisotzkH (120111) suggested a scenario where the 
active fraction of quasars may be a function of BH mass. Cou¬ 
pled with the intrinsic scatter in the M, — cr* relation, this ef¬ 
fect would also produce deviations from the mean M. - cr* 
relation in samples with quasar-based selection. However, 
at low z, there appears to be no significant difference in 
the M, - cr* relation in inactive a nd active galaxies (e.g., 
iWoo et aO2013llBennert et al.l2014l) . The situation, however, 
may be different at high z. 

The “intrinsic” scatter in BH mass around the best-fit re¬ 
gression listed in Table [2] is the combination of the intrinsic 
scatter in the M, - cr* relation and the scatter in SE virial 
BH mass estimates around true BH masses. At face value, 
this suggests that the systematic uncertainty in SE virial BH 
masses is less than ~ 0.4dex, since there must be intrinsic 
scatter in the M. - cr* relation. This statement is incorrect, 
however: because of the luminosity threshold and the way the 
two biases operate, we are not sampling the full distribution 
of objects and thus do not probe the true uncertainty in SE 
virial masses. If we were to select a sample with fixed con¬ 
tinuum luminosity and fixed line width, then the dispersion 
in the sample virial BH masses would be zero, while the un¬ 
certainty in these BH mass estimates would still be ~ 0.5 dex. 
IShenl (120131 Sec 3.3) describes in detail how to distinguish 
these different concepts. In brief, the observed distribution 
in the SE virial masses for a given sample is the convolu¬ 
tion of the underlying BH mass distribution with the lumi¬ 
nosity selection function and the SE mass error distribution. 
Luminosity-threshold quasar samples often show a dispersion 
of virial masse s that can be smaller than the nominal m ass un¬ 
certainty (e.g., IShen et all2008ilShen & KellvBoTl . These 
two are equal only when: (1) the true masses have zero sam¬ 
ple dispersion; and (2) the error in SE virial masses does not 
d epend on luminosity - neither of these are true, as discussed 
in lSheiJd20l3b . 


Similarly, even if we can account for the scatter in SE 
BH masses around true BH masses, the dispersion in true 
BH mass for the luminosity threshold sample is not the in¬ 
trinsic scatter around the mean M. - a, relation, precisely 
due to the Lauer et al. bias: we are sampling the upper en- 
velop e of the M, - cr* relation rather than the full distribu¬ 
tion. iJahnke & Macciol (1201 If) proposed an int e resting sce¬ 
nario, based on an earlier suggestion by iPeng I (120071) . that 
the M. - <7* relation can be naturally produced by the aver¬ 
aging process during mergers of galaxies with initially un¬ 
correlated BH mass and galaxy properties. Such a scenario 
predicts an increasing intrinsic scatter in the M. — cr* rela¬ 
tion with redshift, and is testable with observations. Unfor¬ 
tunately, we cannot test this prediction with our data without 
a more careful forward modeling of the complications just de¬ 
scribed. Nevertheless, the nearly constant scatter in Table [2] 
albeit only measuring the upper envelop of the M. — cr* rela¬ 
tion due to selection effects, seems to favor a non-evolving 
intrinsic scatter in the M m - cr* relation. 

Our obs e rvatio ns are consistent with I Woo et al.l (120061) and 
IWoo et all (120081) : if we look only at their samples at z = 
0.36 and z = 0.57, which are at comparable luminosities (see 
Fig. 121), then there is little evolution. A strong evolution 
was inferred only when they compared to a low-z and low- 
luminosity sample, as shown here as well (Fig. [6]). Given 
the limited dyna mic range in luminos ity in the samples of 
IWoo et all (120061) and I Woo et al.l (H008), this evolution mani¬ 
fested as an offset from the local relation. However, the low- 
Z samples have on average much lower luminosity than the 
high-z samples, so this comparison suffers from the statistical 
biases discussed above. The I Woo et alJ d2006l 120081) samples 
are on average more luminous than our quasars at the rele¬ 
vant redshifts (see Fig. [21, and they do tend to lie in the upper 
envelop of the distribution of our quasars in the Mbh.vu- - <t* 
plane, as shown in Fig. 0] By going to lower luminosities 
at these redshifts, our sample was able to recover some (but 
not all) of the low-mass BHs at fixed cr* that were missed 
in earlier samples. All these observations are consistent with 
the expectation from statistical biases that we demonstrated in 
Fig. [7] with simple models. 

Our interpretation of the observed MBH.vir ~ o'* relation 
at high z is therefore fundamentally different from Woo 
et al., in that we believe the differences seen between high- 
Z and low-z samples are mainly not due to real evolution, 
but largely caused by selection biases and BH mass uncer- 


sized in several recent studies (e.g., Salviander et al. 

2007; 

Lauer et al. 200'|^ Jahnke et al. 

2009 

; Shen & Kelly 

2010; 

Schulze & Wisotzki 2011:; Shen 

2013 

; Salviander & Shields 

2013| Schramm & Silverman 2013; 

Matsuoka et all 201-3; 

Schulze & Wisotzki 2014t Sun et al. 201'®. Indeed, such a 


strong evolution in the slope (or the normalization) of the 
M m - cr* relation to z ~ 1 is inconsistent with the latest cos¬ 
mological hydrody namical simulations dKhandai et alJl2Q14l : 
ISiiacki et alJl20Tl iDeGraf et alJlMl , or independent con¬ 
straints from evolutionary models of the velocity dispersion 
function and accreted BH mass density (IShankar et al.ll20()9b . 
which found no significant evolution in the M. - cr* relation. 


5. CONCLUSIONS 

We have presented a large, homogenous set of new host 
stellar velocity dispersion measurements in a uniformly- 
selected sample of quasars at 0.1 < z < 1, with 46 quasars 
at z > 0.6, a regime rarely explored in previous studies. The 
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sample size and dynamic range in luminosity show signifi¬ 
cant improvement over the best samples studied previously 
(e.g.. IWoo et all[20061120081) . This work was made possible 
by the high S/N coadded spectra of a homogeneous quasar 
sample covering a wide luminosity and redshift range from 
the SDSS-RM project, the first major multi-object reverbera¬ 
tion mapping project to date (IShen et al.ll2015l) . 

We decomposed the spectra into host and quasar compo¬ 
nents using the PCA method, and measured cr* and SE virial 
mass estimates based on the decomposed galaxy and quasar 
spectra, respectively. 

We found the first direct, conclusive evidence that BH mass 
correlates with <r* in quasars at z > 0.3, such that quasar hosts 
with higher cr* have more massive BHs. Similar claims have 
previously been made, but were based on proxies for cr* that 
are much easier to measure in quasar spectra but with signifi- 
cant caveats (such as the width of the [O ill or TO m l line, e.g., 
tSalviander et al.ll2007l : lSalviander & Shieldsii~2013l) . This cor¬ 
relation holds to z ~ 0.8. We do not find significant evolution 
in this relation with our large sample and uniform analysis, up 
to z ~ 1 • 

However, the relation is significantly shallower 

than the local M. — cr* relation (Figs. [3] and |4j. We argued 
(e.g.. Figs. [6] and O that this is most likely due to the sta¬ 
tistical biases induced by the combination of the luminosity- 
threshold quasar selection and the uncertainties in the SE 
virial BH masses. The observed trend of the /W B Rvir — cr* re¬ 
lation with luminosity (both within our sample and combined 
with other samples) appears to confirm this suggestion, but 
clearly more data and tests are needed. It is possible that there 
is still some real evolution to some extent. 

Our work represents a step forward in studies of the evo¬ 
lution of the M, - cr* relation, and demonstrates the potential 
and caveats of using quasars and virial BH mass estimates to 
probe the co-evolution of BHs and galaxies up to high red- 
shift. Our SDSS-RM program has reached a depth and sam¬ 
ple size that are difficult to surpass with traditional programs 
on current 6-10 m telescopes, but future installment of multi¬ 
object RM programs on 6-10 m telescopes or spectroscopy 
with 30 m-class telescopes will be able to provide even fainter, 
and hence less biased, quasar samples to measure host stellar 
velocity dispersion. With future direct RM-based BH masses 
for these quasars (the main goal of the SDSS-RM project), 
we can largely eliminate the luminosity-dependent bias in SE 
virial masses due to the explicit luminosity dependence of 
these mass estimators. It is unlikely that we will obtain a sam¬ 
ple with RM masses and cr* measurements nearly as large as 
the one studied here. Nevertheless, these RM masses will help 
us understand better the systematics in these SE virial masses. 
With fainter quasar samples at high z, we can potentially re¬ 
cover the population of BHs that are undermassive relative to 
their hosts. With these two future advances, we will be able to 
better characterize the relation between BHs and host galax¬ 
ies at high z, which in turn will improve our understanding of 
galaxy formation and BH feedback. 


We thank the anonymous referee for useful comments and 
suggestions, and Scott Tremaine and Andrew Benson for 
their insights. Support for the work of YS was provided by 
NASA through Hubble Fellowship grant number HST-HF- 
51314, awarded by the Space Telescope Science Institute, 
which is operated by the Association of Universities for Re¬ 
search in Astronomy, Inc., for NASA, under contract NAS 5- 


26555. LCH acknowledges support by the Chinese Academy 
of Sciences through grant No. XDB09030102 (Emergence of 
Cosmological Structures) from the Strategic Priority Research 
Program and by the National Natural Science Foundation of 
China through grant No. 11473002. WNB acknowledges sup¬ 
port from NSF grant AST-1108604. BMP acknowledges sup¬ 
port from the NSF through grant AST-1008882. 

Funding for SDSS-III has been provided by the Alfred 
P. Sloan Foundation, the Participating Institutions, the Na¬ 
tional Science Foundation, and the U.S. Department of 
Energy Office of Science. The SDSS-III web site is 
http://www.sdss3.org/. 

SDSS-III is managed by the Astrophysical Research Con¬ 
sortium for the Participating Institutions of the SDSS-III Col¬ 
laboration including the University of Arizona, the Brazilian 
Participation Group, Brookhaven National Laboratory, Uni¬ 
versity of Cambridge, Carnegie Mellon University, University 
of Florida, the French Participation Group, the German Partic¬ 
ipation Group, Harvard University, the Instituto de Astrofisica 
de Canarias, the Michigan State/Notre Dame/JINA Participa¬ 
tion Group, Johns Hopkins University, Lawrence Berkeley 
National Laboratory, Max Planck Institute for Astrophysics, 
Max Planck Institute for Extraterrestrial Physics, New Mex¬ 
ico State University, New York University, Ohio State Univer¬ 
sity, Pennsylvania State University, University of Portsmouth, 
Princeton University, the Spanish Participation Group, Uni¬ 
versity of Tokyo, University of Utah, Vanderbilt University, 
University of Virginia, University of Washington, and Yale 
University. 

REFERENCES 

Adelberger, K. L„ & Steidel, C. C. 2005, ApJ. 627, LI 
Aller, M. C., Richstone, D. O. 2007, ApJ, 665, 120 
Barth, A. J„ Greene, J. E„ & Ho. L. C. 2005, ApJ, 619, L151 
Bennert, V. N„ et al. 2010, ApJ, 708, 1507 

Bennert, V. N„ Auger, M. W„ Treu, T„ Woo, J.-H., & Malkan, M. A. 2011, 
ApJ, 726, 59 

Bennert, V. N., et al. 2014, ApJ, in press, arXiv: 1409.4428 
Benson, A. J. 2010, Physics reports, 495, 33 
Bentz, M. C., et al. 2009, ApJ, 694, L166 

Bernardi, M., Shankar, F., Hyde, J. B., Mei, S., Marulli, F., & Sheth, R. K. 

2010, MNRAS, 404, 2087 
Bolton, A., et al. 2012, AJ, 144, 144 
Bongiorno, A., et al. 2014, MNRAS, 443, 2077 
Cappellari, M„ & Emsellem, E. 2004, PASP, 116, 138 
Canalizo, G., et al. 2012, ApJ, 760, 38 
Cistemas, M., et al. 2011, ApJ, 741, LI 1 

Collin, S., Kawaguchi, T., Peterson, B. M., & Vestergaard, M. 2006, A&A, 
456, 75 

Dasyra. K. M„ et al. 2007. ApJ, 657, 102 
Dawson, K. S., et al. 2013, AJ, 145, 10 
Decarli, R„ et al. 2010, MNRAS, 402, 2453 
DeGraf, C., et al. 2014, arXiv: 1412.4133 
Eisenstein, D. J., et al. 2011. AJ, 142, 72 
Feng, H„ Shen, Y„ & Li, H. 2014, ApJ, 794, 77 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Ferrarese, L., et al. 2001, ApJ, 555, L79 
Gebhardt, K.. et al. 2000a, ApJ, 539, L13 
Gebhardt, K.. et al. 2000b, ApJ, 543, L5 

Graham, A. W., Erwin, P, Caon, N., & Trujillo, I. 2001, ApJ, 563, LI 1 
Graham. A. W., Onken, C. A., Athanassoula, E., Combes, F. 2011, MNRAS, 
412, 2211 

Greene, J. E„ & Ho, L. C. 2006a, ApJ, 641. 117 
Greene, J. E„ & Ho, L. C. 2006b, ApJ, 641, L21 
Grier, C. J., et al. 2013, ApJ, 773, 90 
Gunn, J. E„ et al. 2006, AJ, 131, 2332 
Gultekin, K„ et al. 2009. ApJ, 698, 198 
Haring, N„ & Rix, H.-W. 2004, ApJ, 604, L89 
Hiner, K. D„ et al. 2012, ApJ, 756, 162 
Jahnke, K„ et al. 2009, ApJ, 706, L215 









12 


SHEN ET AL. 


Jahnke, K„ & Maccio, A. V. 2011, ApJ, 734, 92 
Kelly, B. C. 2007, ApJ, 665, 1489 
Khandai, N„ et al. 2014, arXiv: 1402.0888 
Kormendy, J„ & Ho, L. C. 2013, ARA&A, 51, 511 
Lauer, T. R., Tremaine, S.. Richstone, D., & Faber, S. M. 2007, ApJ, 670, 
249 

Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 

Matsuoka, Y., Strauss, M. A., Price, HI, T. N., & DiDonato, M. S. 2014, 

ApJ, 780, 162 

McConnell, N. J„ & Ma, C.-P. 2013, ApJ, 764, 184 
McLeod, K. K„ & Bechtold. J. 2009, ApJ, 704, 415 
Merloni, A., et al. 2010, ApJ, 708, 137 

Moultaka, J., Ilovaisky, S. A., Prugniel, P., & Soubiran, C. 2004, PASP, 116, 
693 

Nelson, C. H„ et al. 2004, ApJ, 615, 652 
Onken, C. A., et al. 2004, ApJ, 615, 645 

Peng, C. Y., Impey, C. D., Ho, L. C., Barton, E. J., & Rix, H.-W. 2006a, ApJ, 
640, 114 

Peng, C. Y., et al. 2006b, ApJ, 649, 616 
Peng, C. Y. 2007, ApJ, 671, 1098 
Peterson, B. M. 2014, Space Science Reviews, 183, 253 
Salviander, S., Shields, G. A., Gebhardt, K., & Bonning, E. W. 2007, ApJ, 
662, 131 

Salviander, S„ & Shields, G. A. 2013, ApJ, 764, 80 
Schramm, M., & Silverman, J. D. 2013, ApJ, 767, 13 
Schulze, A., & Wisotzki, L. 2011, A&A, 535, 87 
Schulze, A., & Wisotzki, L. 2014, MNRAS, 438, 3422 
Shankar, F„ Bemardi, M., & Haiman, Z. 2009, ApJ, 694, 867 
Shen, J., Vanden Berk, D. E„ Schneider, D. P„ & Hall, P. B. 2008b, AJ, 135, 
928 


Shen, Y. 2013. BASI, 41, 61 
Shen, Y„ et al. 2011, ApJS, 194, 45 
Shen, Y„ et al. 2015, ApJS, 216, 4 

Shen, Y„ Greene, J. E„ Strauss, M. A., Richards, G. T„ & Schneider, D. P. 
2008a, ApJ, 680, 169 

Shen, Y„ & Ho, L. C. 2014, Nature, 513, 210 

Shen, Y„ & Kelly, B. C. 2010, ApJ, 713, 41 

Shen, Y„ & Kelly, B. C. 2012, ApJ, 746, 169 

Sijacki, D„ et al. 2014, MNRAS, arXiv: 1408.6842 

Smee, S. A., et al. 2013, AJ, 146, 32 

Sun, M., et al. 2014, ApJ, in press, arXiv: 1502.01025 

Targett, T. A., Dunlop, J. S., & McLure, R. J. 2012, MNRAS, 420, 3621 

Tremaine, S., et al. 2002, ApJ, 574, 740 

Treu, T„ Malkan, M. A., & Blandford, R. D. 2004, ApJ, 615, L97 
Treu, T„ Woo, J.-H., Malkan, M. A., & Blandford, R. D. 2007, ApJ, 667, 117 
Valdes, F., Gupta, R.. Rose, J. A., Singh, H. P., & Bell, D. J. 2004, ApJS, 
152,251 

Vanden Berk, D. E„ et al. 2006, AJ, 131, 84 
Vestergaard, M., & Peterson, B. M. 2006, ApJ, 641, 689 
Wang^J., et al. 2009, ApJ, 707, 1334 
Watson, L. C., et al. 2008, ApJ, 682, L21 
Woo, J.-H., et al. 2013, ApJ, 772, 49 

Woo, J.-H., Treu, T„ Malkan, M. A., & Blandford, R. D. 2006, ApJ, 645, 900 
Woo, J.-H., Treu, T„ Malkan, M. A., & Blandford, R. D. 2008, ApJ, 681, 925 
Yip, C. W„ et al. 2004a, AJ, 128, 585 
Yip, C. W„ et al. 2004b, AJ, 128, 2603 



